Thermoelectric and thermal rectification properties of quantum dot junctions 
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The electrical conductance, thermal conductance, thermal power and figure of merit (ZT) of 
semiconductor quantum dots (QDs) embedded into an insulator matrix connected with metallic 
electrodes are theoretically investigated in the Coulomb blockade regime. The multilevel Anderson 
model is used to simulate the multiple QDs junction system. The charge and heat currents in the 
sequential tunneling process are calculated by the Keldysh Green function technique. In the linear 
response regime the ZT values are still very impressive in the small tunneling rates case, although 
the effect of electron Coulomb interaction on ZT is significant. In the nonlinear response regime, 
we have demonstrated that the thermal rectification behavior can be observed for the coupled QDs 
system, where the very strong asymmetrical coupling between the dots and electrodes , large energy 
level separation between dots and strong interdot Coulomb interactions are required. 



I. INTRODUCTION 

Due to energy and environment issues, it becomes im- 
portant to understand the thermal properties of mate- 
rials. Recently many efforts are to seek efficient ther- 
moelectric materials because there exist potential ap- 
plications of solid state thermal devices. 1-9 Neverthe- 
less, the optimization of thermoelectric properties of ma- 
terials is extremely difficult, since the figure of merit 
(ZT = S 2 G e T/n) depends on Seebeck coefficient (S), 
electrical conductance (G e ) and thermal conductance (n) 
of the material. Tuning one of these physical quantities 
will unavoidably alter the other because they are closely 
related. 7 

Several methods were proposed to realize the en- 
hancement of ZT, 2 one of them is to reduce the sys- 
tem dimensionality. 8 Bi2Te3/Sb2Te3 superlattices 3 , sil- 
icon quantum wires 4 and PbSeTe based quantum dot 
(QD) superlattices 5 were experimentally demonstrated 
to show much higher ZT values when compared with 
their corresponding bulk materials. A zero-dimension 
QD system was predicted to have more pronounced en- 
hancement in thermoelectric efficiency due the reduced 
dimensionality 9 . Experimentally, it has been shown 5 
that the performance of PbSeTe QDs can reach a very im- 
pressive ZT value of 2. Nevertheless, a ZT value higher 
than 3 has never been reported. Note that the high- 
est ZT value is near 1 for conventional bulk materials. 1 
Systems with ZT value larger than 3 may find appli- 
cation in making home refrigerators, replacing the exist- 
ing compressor-based refrigerators. In addition, they can 
used in electrical power generators. 2 

In order to seek a large ZT value, a single molecular 
QD weakly linked to electrodes was proposed to exhibit 
an extremely large ZT value in the Coulomb blockade 
regime. 10 However, reference [10] did not take into ac- 
count the molecular vibrations. For a molecular junction, 
the coupling strengthes between localized electrons and 
vibration modes are very strong. 11-17 Due to multiple 
phonon assisted processes arising from strong electron 



phonon interactions, it is expected that ZT values will 
be suppressed by molecular vibrations. Apart from that, 
such a molecular junction is difficult to integrate with 
current silicon based electronics. Therefore, we propose 
to use a thermoelectric device made of semiconductor 
QDs embedded into amorphous insulator which has low 
heat conductivity. The studied system is shown in Fig. 
1. In addition, a nanoscale vacuum layer is inserted to 
block the heat current delivered by phonon carriers, al- 
though it would be a challenging task to keep the vacuum 
layer thin enough to allow sufficient electron tunneling. 
The vacuum layer considered here can be realized using 
the technology similar to that used in liquid crystal dis- 
play implementation where a vacuum layer is inserted for 
blocking the heat generated by the light source. 

The key applications of thermoelectric devices include 
solid state refrigerators and electrical generators. In solid 
state refrigerators (electrical generators), one needs to 
remove (generate) large amount of heat current (charge 
current). Consequently, a high QD density is required 
for realistic applications. A single level Anderson model 
can be used to simulate such a system adequately in the 
dilute QD-density limit. 10 However, for the high QD- 
density system, one needs to consider the effect of in- 
terdot Coulomb interactions and electron hopping effect. 
When QDs are embedded in an insulator matrix having a 
high potential barrier, electron hopping among dots can 
be neglected. However, it is hard to avoid the interdot 
Coulomb interactions due to its long-range tail. 

In this paper, we investigate the effect of interdot 
Coulomb interactions on the thermoelectric properties in 
the linear and nonlinear response regimes via a multi- 
level Anderson model 18,19 . We found that the inter- 
dot Coulomb interactions would suppress the ZT values 
and play a crucial role in determining the thermal rec- 
tification behavior. The electrical conductivity, thermal 
power, thermal conductivity and figure of merit were typ- 
ically calculated in the linear response regime, while cru- 
cial applications of thermal devices in thermal rectifiers 
and transistors require the understanding of the thermo- 



2 



electric properties in the nonlinear response regime. 20 ' 21 
The thermal rectifiers can be used in solar energy stor- 
age and many other applications. Therefore, it is im- 
portant to take into account the thermoelectric effects in 
the nonlinear regime. Here, we demonstrate that coupled 
QDs can exhibit pronounced thermal rectification behav- 
ior. Although the mechanism of thermal rectification for 
QD junctions is similar to that of charge current, the 
heat current is generated by temperature gradient and 
the consequent electrochemical potential. It is the non- 
linear relation between the applied temperature gradient 
and the electrochemical potential that leads to enhance 
thermal rectification behavior. 



II. FORMALISM 

A schematic diagram of the system of concern is shown 
in Fig. 1. The Hamiltonian of the system can be de- 
scribed by a multi-level Anderson model: 
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where a\ a ^ {ak.a,p) creates (destroys) an electron of mo- 
mentum k and spin a with energy in the (3 metallic 
electrode. d\ a (c^ i<7 ) creates (destroys) an electron with 
the ground-state energy Eg in the £th QD , Ug and Ugj 
describe the intradot Coulomb interactions and the inter- 
dot Coulomb interactions, respectively. Vk,p,g describes 
the coupling between the band states of electrodes and 
the QD levels. We have ignored the excited levels of 
QDs, assuming that the energy level separation between 
the ground state and the first excited state within each 
QD is much larger than intradot Coulomb interactions 
Ug and thermal energy k R T, where T is the temperature 
of concern. We have also ignored the inderdot hopping 
terms due to the high potential barrier separating QDs. 
The key effects included are the intradot and interdot 
Coulomb interactions and the coupling between the QDS 
with the metallic leads. 

Using the Keldysh-Green's function technique, 22 ' 23 the 
charge and heat currents leaving electrodes can be ex- 
pressed as 



where the transmission factor is 7/fe) = w^ttt4^t-\ ■ 

./Li?(e) = /z,(e) - /fl(e), where f L {R){t) = 
l/[ e (e-»L(R))/k B T HR) + X j ig the Fermi distribution 
functions for the left (right) electrode . The chemical 
potential difference between these two electrodes is 
related to the bias difference via he — = eAV. Ep 
is the Fermi energy of electrodes. Tg^{e) and Tg tR (e) 
= 27r^ k |V^,/3,k| 2 S(e — ek)] denote the tunneling 
rates from the QDs to the left and right electrodes, 
respectively, e and h denote the electron charge and 
Plank's constant, respectively. For simplicity, these 
tunneling rates will be assumed energy- and bias- 
independent. Therefore, the calculation of tunneling 
current and heat current is entirely determined by 
the spectral function, A(e) = ImG r lcr {e), which is the 
imaginary part of the retarded Green's function G\ a (e) . 
The expression of retarded Green function is given 
by[18,19] 
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where n denotes the number of coupled QDs in each cell 
considered. H m denotes the sum of Coulomb interactions 
seen by a particle in dot i due to other particles in the 
dot j(j ^ £), which can be occupied by zero, one or two 
particles. p m denotes the probability of such configura- 
tions. For a three-QD cell (£ ^ j ^ j'), there are nine 
(3 x 3) configurations, and the probability factors become 

Pl = djCLf, p 2 = bjCLy, p 3 = djbji, Pi = Cjdj', p 5 = Cj'dj, 

Pq = bjbf, pi = cjbj>, ps = Cfbj, and p 9 = CjCj>, where 
Oj = 1 - (Nj,a + Nj,_ a ) + Cj, bj = {Nj.a + N jt _ a ) - 2c„ 
and cj — (riji-o-nj^) is the intradot two-particle corre- 
lation function. Nj, a is one particle occupation num- 
ber. Interdot Coulomb interaction factors are III = 0, 

n 2 = u t ,i, n 3 = u Lr , n 4 = 2U t}j , n 5 = 2u e . r , 

n 6 = U e ,j + U e ,j>, n 7 = 2U e ,j + U e , r , n 8 = 2U t . r + U ej , 
and n 9 = 2U e ,j + 2U e ,j>. Tg = (T tiL + T LR )/2 arises 
from the self-energy due to the weak coupling between 
the QDs with metallic leads, where the real part of self 
energy is ignored. Such a self energy (ignoring the ef- 
fect of electron Coulomb interactions) is adequate within 
the Coulomb blockade regime, but it does not capture 
the Kondo effect. The sum of probability factors p m for 
all configurations is equal to 1, reflecting the fact that 
G\ CT (e) satisfies the sum rule. 

According to the expression of retarded Green's func- 
tion of Eq. (4), we need to know the single-particle 
and two-particle occupation numbers, Ng tlJ (Ng_ a ) and 
Ne,e — eg, which can be obtained by solving the follow- 
ing equations self-consistently. 
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The values of N^ a and ct are restricted between and 1. 
The expression of two particle retarded Green function 
of Eq. (6) is 
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III. LINEAR REGIME 

In the linear response regime, Eqs. (2) and (3) can be 
rewritten as 



. AV AT 

Je - i-ll-y- + 
y - 1-21— + 



(7) 



where AT = Tl — T# is the temperature difference across 
the junction. Coefficients in Eq. (7) are given by 
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Here T(e) = - I^TJiFTO ^(OlAy=o,Ar=o 
and /(e) = i/[ e (e-E F )/k B T + ^ Note that thc Qn _ 
sager relation £12 = £21 is preserved. Based on Eq. 
(7), thc charge current can be generated by the volt- 
age difference and temperature gradient. If the system 
is in an open circuit, the electrochemical potential will 
form in response to a temperature gradient; this electro- 
chemical potential is known as the Seebeck voltage (See- 
beck effect). Seebeck coefficient ( the amount of voltage 
generated per unit temperature gradient) is defined as 
S = = — T§77- I R terms of the Seebeck coefficient, 
the electron thermal conductance is n e — (^ — veils' 2 ). 
To judge whether the system is able to generate power 
or refrigerate efficiently, we need to evaluate the figure 
of merit, ZT = S 2 G e T/k, where G e = is the elec- 

trical conductance and k = n e + K p h is the thermal con- 
ductance. K p h denotes the thermal conductance due to 



the phonon contribution. . For a system with an effi- 
cient thermoelectric properties we want ZT as high as 
possible. This implies that we desire a system with high 
Seebeck coefficient, high electrical conductance and low 
thermal conductance. The thermal conductance arising 
from phonons can be neglected (k — n e ) in our proposed 
system because the vacuum layer can block the heat cur- 
rent carried by phonons effectively. 

Although Eq. (4) can be employed to calculate the 
charge current and heat current of a junction system with 
arbitrary QD number, 18 ' 19 here we use the three-QD ex- 
ample to investigate the effect of interdot Coulomb inter- 
action on the figure of merit, ZT. As mentioned above, 
ZT depends on the electrical conductance G e , Seebeck co- 
efficient S and electron thermal conductance n e . There- 
fore, it is difficult to calculate the exact solution of ZT 
for arbitrary parameters. For simplicity, we have ignored 
the QD size fluctuations and assumed all QDs have the 
same ground-state energy, Eg = E g in the evaluation 
of ZT. Thc QD size fluctuations will become important 
in the consideration of heat current rectification below. 
TClosed form expressions for the coefficients defined in 
Eqs. (8)-(ll) exist within the small tunneling-rate limit 
(i.e. ( e _ E T y? +r i! ji can be approximated by ir5(e — E g )) 
and no electron Coulomb interaction. We obtain 

£11 =a /cosh 2 (A/(2fc B T)), 
£12 = C21 = ai/cosh 2 (A/(2/c B T)), 
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ductance K e = (^# — AiS" 2 ) vanishes, whereas the elec- 
trical conductance of G e — y£n and the thermal power 
of S = — remain finite. This indicates that system 
ZT diverges as T approaches zero. This is the so called 
"Carnot efficiency". 10 

Closed-form expressions for these coefficients for finite 
r in the non-interacting case have also been derived in 
terms of trigamma functions. 10 However, the complicated 
trigamma functions do not simplify the expression of ZT 
and make it difficult to elucidate mechanisms for optimiz- 
ing ZT. Therefore, we numerically calculate the figure of 
merit ZT with and without the Coulomb interactions. 
We first consider the case of symmetrical tunneling rates 
{Tl = Yr = lmeV). The consideration of asymmet- 
rical tunneling rates is not important for the linear re- 
sponse regime, but it is crucial for the nonlinear response 
regime. 20 ' 21 which we shall address in the next section. 
ZT as a function of temperature for various values of A in 
the absence of interdot Coulomb interactions is shown in 
Fig. 2. Solid lines and dotted lines denote cases without 
and with intradot Coulomb interaction (U = 125r), re- 
spectively. Note that all energies are measured in terms 
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of r through out this article. We see that the solid lines 
merge with the dashed lines at low temperatures. This 
indicates that the effect of intradot Coulomb interaction 
on ZT can be ignored when U ksT. Such a result can 
be understood as follows. When interdot Coulomb inter- 
actions Utj = Uds vanish, the retarded Green function 
consists of two branches 

The second branch has a negligible contribution due to 
the vanishing factor e xp~( Eg+u ~ EF ^ kBT (when U 3> 
ksT and A > 0) which appears in Eqs. (8)-(ll). The 
factor (1—Ni >a ) in the first branch of Eq. (12) only affects 
the coefficients (£n, £12, £21 and £22), but not their ra- 
tios. This explains why ZT(U = 0) « ZT(U > k B T). 
The reduction of ZT at finite U (for instance U/ksT = 5) 
can be understood as follows. In the small tunneling 
rate limit, we find that S 2 G e oc T and n e oc T for finite 
U. This is different from the behavior, S 2 G e oc T and 
n e oc T 2 in the absence of U. Consequently, the reduc- 
tion of ZT is observed in Fig. 2. 

For a thermal electric device with high QD density, 
the interdot Coulomb interactions are also important. 
Fig. 3 shows the ZT value as a function of tempera- 
ture for a three-QD cell for various QD configurations 
with A = 30 T and U = 125 T. Dotted line de- 
notes the case of dilute QD density. As a result of a 
large separation between QDs, the interdot Coulomb in- 
teractions are negligible (Ue,j = 0). Dashed lines de- 
notes the case where dot A and dot B are close to 
each other (Uab = 45r), but dot C is far from them 
(Uac — Ubc = 0). Dot-dashed line denotes the case 
with Uab = U A c = 45r but U B c = 0. Solid line de- 
notes the case with Uab = Uac = Ubc = 45r. The 
results of Fig. 3 indicate considerable reduction of ZT 
at high temperatures due to the interdot Coulomb in- 
teractions (proximity effect). However, the proximity 
effect on ZT can be ignored when Uej/k B T >• 1. In 
general, the maximum values of interdot Coulomb in- 
teractions are one-half of intradot Coulomb interactions. 
For silicon QDs embedded in S1O2, the Si QDs with 
5 — 10 nm diameters have the intradot Coulomb inter- 
action strengthes between 100 — 150 meV. Therefore, the 
condition of U(j/k B T ^S> 1 is not easy to be satisfied at 
room temperature. 

According to the results of Fig. 2, the system ZT can 
be tuned by the A value. In Fig. 4, we plot ZT as a 
function of A for the three-QD system with and without 
Coulomb interactions for various temperatures. For the 
nonintcracting case (thin solid line with mark), the max- 
imum ZT value occurred at A max = 2Ak B T, which has 
been pointed out in Ref. 10. However, for the case with 
finite electron Coulomb interactions, we found A maXt i = 
3.1k B T, A max , 2 = 3k B T and A max , 3 = 2.8k B T. It is 
worth noting that the maximum ZT values for different 
temperatures still reach ZT max > 3, which are very en- 



couraging values. However, we have not considered the 
QD size fluctuations, defects between metallic electrodes 
and insulators, and clcctron-phonon interactions. In or- 
der to include these affects fully, we phenomenologically 
replace the imaginary part of retarded Green function of 
Eq. (4) by Tg + Ti e . This means the total level-width 
is expressed as the sum of elastic and inelastic widths. 
Fig. 5 shows the inelastic scattering effect on ZT. For 
Ti e = 3r, the maximum ZT value becomes smaller than 
4. The results of Fig. 5 indicate that the suppression of 
ZT max resulting from the inelastic scattering is serious. 

To further understand the results of Figs. 3 and 4, 
we analyze the electron conductance G e , thermal power 
S and electron thermal conductance n e of the system. 
Fig. 6 shows G e , S and n e as functions of temperature 
for a thrcc-QD cell for various configurations: Uab = 
Uac = Ubc = 45r (dotted curves), U A b = U A c = 45r 
and Ubc = (dashed curves), and Uab = 45r and 
Uac — Ubc = (solid curves), which correspond to 
strong, medium, and weak proximity effects. We noticed 
from Fig. 6(b) that the thermal power (S) is not sensitive 
to the proximity effect, which means the proximity effects 
on £12 and C\\ are similar. Thus, the ZT behavior at 
high temperature shown in Fig. 3 is mainly attributed 
to K e and G e . n e is enhanced [see Fig. 6(c)], but G e 
is suppressed [see Fig. 6(a)] at high temperature when 
the proximity effect increases. This explains why ZT is 
suppressed at high temperature with increasing proxim- 
ity effect. The maximum absolute value of S appears at 
near ksT = 6T, whereas the maximum ZT value shown 
in Fig. 3 appears between ksT = 10T and fc^T = 15r. 
So the temperature dependence of ZT is similar to that 
of the electrical conductivity G e , meaning that S 2 T '/ ' n e 
has a weak temperature dependence. 

In Fig. 4, we have tuned A from to 120r. This 
implies that the energy levels of QDs are shifted away 
from the Fermi energy of electrodes. ZT becomes small 
when A/ksT 3> 1. We can apply a gate voltage (V g ) to 
move Eg relative to Ep. In. Fig. 7 we plot G e , S and n e 
as functions of the gate voltage (V^) for various tempera- 
tures with A fixed at 30r. The electric conductance (G e ) 
clearly exhibits a Coulomb oscillation arising from the in- 
tradot and interdot Coulomb interactions. The first three 
peaks of G e result from the resonant channels of poles at 
Eg, Eg + Uds, and E g + 2Ud s for e. Other peaks can 
be readily identified by the resonant channels of retarded 
Green function of Eq. (4). Note that the resonant chan- 
nel of e = Eg + U (= 155r) is seriously suppressed due 
to the fact that all three QDs are filled with one electron 
at that gate voltage. The Coulomb oscillatory behavior 
of G e becomes smeared at higher temperatures. 

The thermal power (S) exhibits a sawtooth-like shape 
with respect to gate voltage, which is consistent with 
the experimental observation. 24 The sawtooth-like shape 
was also theoretically reported in the metallic single elec- 
tron transistor, 25 ' 26 where the charging energies arc ho- 
mogeneous. In Refs. 25 and 26, a model based on the 
rate equations was adopted. The thermal power S can 
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be tuned from negative to positive values. When the 
Fermi energy matches a resonant channel (G e reaches 
a maximum value) the thermal power vanishes, since 
£12 = 0. When the Fermi energy of electrodes is in 
the middle of two resonant channels, S also vanishes. 
Zero thermal power indicates that the current arising 
from temperature gradient can be self-consistently bal- 
anced without electrochemical potential. The behavior 
of n e is much more complicated than that of G e , since 
K e = {^pT — £-nS 2 ) which consists of C\\, £12, and 
£ 2 2 ■ Since n e is positive definite, we obtain the relation 
%r > £11 S 2 . When thermal power vanishes, n e — 
Based on the results shown in Fig. 7, the optimized 
ZT value does not match either the maximum G e (good 
conductor) or the minimum G e (poor conductor). The 
largest value for ZT is obtained midway between the good 
and poor conductors as illustrated in Fig. 8 for k B T = 2T 
and A = 30r. So far, our discussion is limited to the 
linear response regime with AT/T <C 1. Some function- 
alities of thermal electric devices require that the applied 
temperature bias AT violates the AT/T <C 1 condition. 
In the following study, the thermoelectric properties of 
QD junctions are investigated in the nonlinear response 
regime. 



IV. NONLINEAR REGIME 

Schcibncr and coworkers experimentally reported the 
thermal power of the two-dimensional electron gas in 
QD under high magnetic fields in the linear response 
regime. 27 Few theoretical works have reported the ther- 
mal properties of QD junctions in the nonlinear response 
regime. 28 Ref. 28 theoretically studied the thermal power 
in the Kondo regime based on one-level Anderson model. 
Here, we study the thermal electric effect of multiple QD 
junction in the Coulomb blockade effect in the nonlinear 
regime. We show that in the nonlinear regime, the ther- 
mal rectification behavior can become quite pronounced. 
Records of thermal rectification date back to 1935 when 
Starr discovered that copper oxide/copper junctions can 
display a thermal diode behavior. 29 Recently, thermal 
rectification effects have been predicted to occur in one 
dimensional phonon junction systems. 30-34 

To study the direction-dependent heat current, we 
let T L = T + AT/2 and T R = T - AT/2, where 
Tq = (Tl + T R )/2 is the equilibrium temperature of two 
side electrodes and AT = T L — T R is the temperature 
difference. Because the electrochemical potential differ- 
ence, eAV yielded by the thermal gradient could be sig- 
nificant, it is important to keep track the shift of the 
energy level of each dot according to eg = Ee + r/iAV/2, 
where rji is the ratio of the distance between dot I and 
the mid plane of the QD junction to the junction width. 
Here we set r] B = Vc = 0. A functional thermal rectifier 
requires a good thermal conductor for AT > 0, but a 
poor thermal conductor for AT < 0. Based on Eqs. (2) 
and (3), the asymmetrical behavior of heat current with 



respect to AT requires not only highly asymmetric cou- 
pling strengthes between the QDs and the electrodes but 
also strong electron Coulomb interactions between dots. 
To investigate the thermal rectification behavior, we have 
numerically solved Eqs. (2) and (3) for multiple-QD junc- 
tions involving two QDs and three QDs for various system 
parameters. We first determine AV by solving Eq. (2) 
with J e = (the open circuit condition) for a given AT, 
Tq and an initial guess of the average one-particle and 
two-particle occupancy numbers, Ng and ci for each QD. 
Those numbers are then updated according to Eqs. (5) 
and (6) until self-consistency is established. Once AV is 
solved, we then use Eq. (3) to compute the heat current. 

Fig. 9 shows the heat currents, occupation numbers, 
and differential thermal conductance (DTC) for the two- 
QD case, in which the energy levels of dot A and dot B 
are Ea = Ep — AE/5 and Eb = Ep + a B AE, where as 
is tuned between and 1. We have adopted AE = 200r, 
which is used to describe the energy level fluctuation 
of QDs. The heat currents are expressed in units of 
Qo = T 2 /(2h) through out this article. The intradot 
and interdot Coulomb interactions used are Ue — 30k B T 
and Uab — mk B To. The tunneling rates are Tar = 0, 
r al = 2T, and Tbr = T B l = T. k B To is chosen to be 
25r throughout this article. Here, T = (T AL + T AR )/2 
is the average tunneling rate in energy units, whose typ- 
ical values of interest are between 0.1 and 0.5 meV. The 
dashed curves are obtained by using a simplified expres- 
sion of Eq. (3) in which we set the average two particle 
occupation in dots A and B to zero (resulting from the 
large intradot Coulomb interactions) and taking the limit 
that r <C k B T$ so the Lorentzian function of resonant 
channels can be replaced by a delta function. We have 

Q/ 7B = ir(l-N B )[(l-2N A )(E B -E F )f LR (E B ) 

+ 2N A (E B + Uab - E F )f LR {E B + U AB )](13) 

Here N A ( B ) is the average occupancy in dot A(B). There- 
fore, it is expected that the curve corresponding to 
Eb = Ep + 4AE/5 obtained with this delta function 
approximation is in good agreement with the full solu- 
tion, since E B is far away from the Fermi energy level. 
For cases when E B is close to Ep, the approximation is 
not as good, but it still gives qualitatively correct behav- 
ior. Thus, it is convenient to use this simple expression to 
illustrate the thermal rectification behavior. The asym- 
metrical behavior of Na with respect to AT is mainly 
resulted from the condition Tar = and T A p = 2T. 
The heat current is contributed from the resonant chan- 
nel with e = E B , because the resonant channel with 
e = E B + U AB is too high in energy compared with Ep. 
The sign of Q is determined by /lr{Eb), which indi- 
rectly depends on Coulomb interactions, tunneling rate 
ratio and QD energy levels. The rectification behavior 
of Q is dominated by the factor 1 — 2N A , which explains 
why the energy level of dot-A should be chosen below 
Ep and the presence of interdot Coulomb interactions is 
crucial. The negative sign of Q in the regime of AT < 
indicates that the heat current is from the right elec- 
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trode to the left electrode. We define the rectification 
efficiency as tjq = (Q(+AT) - \Q(—AT)\)/Q(+AT). We 
obtain r] Q (AT = 30T) = 0.86 for E B = E F + 2AE/5 and 
Vq(AT = 30r) = 0.88 for E B = E F + AAE/5. Fig. 9(c) 
shows DTC in units of Qoks/T. It is found that the rec- 
tification behavior is not very sensitive to the variation of 
Eb- DTC is roughly linearly proportional to AT in the 
range — 20r < fc^AT < 20r. In addition, we also find a 
small negative differential thermal conductance (NDTC) 
for E B = E F + 4AE/5. Similar behavior was reported 
in the phonon junction system. 35 

Note that the mechanism of thermal rectification is 
similar to the charge current rectification. However, the 
heat current is yielded by the temperature bias and the 
electrochemical potential. In particular, the electrochem- 
ical potential is a highly nonlinear function of the temper- 
ature bias, which has never been reported for quantum 
dot junctions. Consequently, it is not straightforward to 
reveal the behavior of heat current with respect to the 
temperature bias. The manifested difference between the 
heat current and the charge current is that the origin of 
NDTC is different from that of negative differential con- 
ductance (NDC). The NDC of charge current requires the 
upper energy levels with the shell-filling condition, which 
was discussed in Refs. 18 and 19. For NDTC, it only ap- 
pears in the lower level with shell-filling condition. Fig. 
10 shows the rectification efficiency as a function AT for 
two different values of Eb- The rectification efficiency 
vanishes when UbAT/T <C 1. This implies that it is 
difficult to judge the rectification effect in the linear re- 
sponse regime of AT/T <C 1. Although the two-dot case 
can reach a high rectification efficiency, the heat current 
should be enhanced from the application point of view. 

Fig. 11 shows the heat current, differential thermal 
conductance and thermal power as functions of tem- 
perature difference AT for a three-QD case for various 
values of T A r, while keeping T B (c),R = r s(C),i = r - 
Here, we adopt i]a = \Tal — Tar\/(2T) instead of fix- 
ing t\a at 0.3 to reflect the correlation of dot position 
with the asymmetric tunneling rates. We assume that 
the three QDs are roughly aligned with dot A in the 
middle. The energy levels of dots A, B and C are cho- 
sen to be Ea — Ep — AE/5, E B = Ep + 2AE/5 
and E c = E F + 3AE/5. U AC = Uba = 15fc B T , 
Ubc = 8fc B T , Uc = 30k B T 0} and all other parameters 
are kept the same as in the two-dot case. The thermal 
rectification effect is most pronounced when Tar = 0. 
as seen in Fig. 12(a). (Note that the heat current is 
not very sensitive to Use)- In this case, we obtain a 
small heat current Q = 0.068Q at AT = -30r, but a 
large heat current Q = 0.33Qo at AT = 30r and the 
rectification efficiency t]q is 0.79. However, the heat cur- 
rent for Tar, = is small. For Tar = O.ir, we obtain 
Q = 1.69Q at AT = -30r,Q = 5.69Q at AT = 30r, 
and T]q = 0.69. We see that the heat current is sup- 
pressed for AT < with decreasing Tar- This im- 
plies that it is important to blockade the heat current 
through dot A to observe the rectification effect. Very 



clear NDTC is observed in Fig. 11(b) for the Far = 0.1T 
case, while DTC is symmetric with respect to AT for the 
Far = T A l case. 

From the experimental point of view, it is easier to 
measure the thermal power than the direction-dependent 
heat current. The thermal power as a function of AT is 
shown in Fig. 11(c). All curves except the dash-dotted 
line (which is for the symmetrical tunneling case) show 
highly asymmetrical behavior with respect to AT, yet it 
is not easy at all to judge the efficiency of the rectification 
effect from S for small |AT| (k B \AT\/T < 10). Thus, it 
is not sufficient to determine whether a single QD can act 
as an efficient thermal rectifier based on results obtained 
in the linear response regime of AT /To <C l. 27 According 
to the thermal power values, the electrochemical poten- 
tial eAV can be very large. Consequently, the shift of 
QD energy levels caused by AV is quite important. To 
illustrate the importance of this effect, we plot in Fig. 12 
the heat current for various values of Ec for the case with 
Tar = 0, Ubc = lOfc^To and tja = 0.3. Other param- 
eters are kept the same as those for Fig. 11. The solid 
(dashed) curves are obtained by including (excluding) the 
energy shift t]aAV/2. It is seen that the shift of QD en- 
ergy levels due to AV can lead to significant change in the 
heat current. It is found that NDTC is accompanied with 
low heat current for the case of Ec = Ep+AE/5 [see Fig. 
12(b)]. Even though the heat current exhibits rectifica- 
tion effect for E c = E F + AE/5 and E c = E F + 3A£/5, 
the thermal power has a very different behavior. From 
Figs. 11(c) and 12(c), we see that the heat current is 
a highly nonlinear function of electrochemical potential, 
AV. Consequently, the rectification effect is not straight- 
forwardly related to the thermal power in this system. 

Because the position distribution fluctuation is com- 
mon for QDs, we investigate the interdot Coulomb inter- 
actions on the rectification effect. Fig. 13 shows the heat 
current, electrochemical potential and occupation num- 
ber as functions of AT for various values of U ac with 
Ec = Ep + AE/5. Other parameters are the same as 
in Fig. 12. When Uac — 0, the rectification efficiency 
is suppressed seriously. The residue rectification mainly 
arises from the correlation between dot A and dot B. 
Such results indicate that it is crucial to control the QD 
position in the implementation of QD thermal rectifiers. 
We find that the electrochemical potential is not signif- 
icantly changed when Uac decreases, whereas the heat 
current has a considerable variation. Fig. 13(c) shows 
the occupation numbers of dots A and C. Nb are ignored 
due to their energy levels being far away from the Fermi 
energy level. It is expected that Na is not sensitive to 
the decrease of Uac- increases so much when Uac 
decreases since the main resonant channels of dot C are 
dominated by Ec rather than the combination of Ec and 
Ec + Uac- Fig. 13(c) reveals that the serious suppression 
of rectification efficiency of dot-dashed line shown in Fig. 
13(a) is mainly attributed to the heat current through 
dot C. We once again investigate the rectification effi- 
ciency for three-dot case. Fig. 14 shows the rectification 
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efficiency as function of AT. Ali other parameters are 
the same as those of Fig. 12. The rectification efficiency 
increases with increasing temperature bias. However, t)q 
is not sensitive to the energy level of dot C. 

Comparing the heat current of the three-dot case 
(shown in Figs. 11 and 12) to the two-dot case (shown in 
Fig. 9), we find that the rectification efficiency is about 
the same for both cases (shown in Figs. 10 and 14), while 
the magnitude of the heat current can be significantly en- 
hanced in the three-dot case. For practical applications, 
we need to estimate the magnitude of the heat current 
density and DTC of the IQV junction device in order 
to see if the effect is significant. We envision a thermal 
rectification device made of an array of multiple QDs 
(e.g. thrce-QD cells) with a 2D density N 2c i = 10 n cm~ 2 . 
For this device, the heat current density versus AT is 
given by Figs. 11 and 12 with the units Qo replaced 
by N2dQo, which is approximately 965W/m 2 if we as- 
sume T = 0.5meV. Similarly, the units for DTC becomes 
^dfcsQo/r, which is approximately 34W / /m 20 K. Since 
the phonon contribution can be blocked by the vacuum 
layer in our design, this device could have practical appli- 
cations near 140°K with (fcsTo « 12.5meV). If we choose 
a higher tunneling rate T > lmeV and Coulomb energy 
> 300meY (possible for QDs with diameter less than 1 
nm) , then it is possible to achieve room-temperature op- 
eration. It is worth pointing out that if the vacuum layer 
is replaced by a typical phonon glass, such as SiC>2, which 
has a thermal conductivity of K p h = 1.5W/m°K [36] at 
room temperature, the heat current carried by phonons 
across a 10 nm junction with a temperature bias of 1°K 
would be around 1.5 x 10 8 W/m 2 . This would completely 
dominate over the thermal electric effect considered here 
(by six orders of magnitude). Therefore, unless a vacuum 
layer is inserted, the term K p h will play a dominant role. 



V. SUMMARY AND CONCLUSIONS 

We have theoretically investigated the effect of intradot 
and interdot Coulomb interactions on the figure of merit 
(ZT) and thermal power (S) of multiple QD junction sys- 
tem in the sequential tunneling process. The ZT val- 
ues at high temperatures are significantly suppressed by 
the intradot as well as interdot Coulomb interactions. 
The optimization of ZT depends not only on tempera- 
ture but also on the detuning energy (A = E g — Ep). 
It is worth noting that inelastic scattering effect arising 
from QD size fluctuations, defects and electron-phonon 
interactions will lead to considerable reduction to the ZT 
values. Electrical conductance and thermal power exhibit 
Coulomb oscillatory behavior and the sawtooth-like be- 
havior with respect to the gate voltage. The largest value 
for ZT is obtained midway between good and poor con- 
ductors. Apart from the results of linear response, the 
heat rectification effect can be observed for multiple QD 
junctions in the nonlinear response regime. In contrast 
to the heat rectification of phonon junction system, the 
heat current is carried by electrons in the multiple QD 
junction system and large electrochemical potentials can 
be established by the temperature gradient to generate 
electrical power. 
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Figure Captions 

Fig. 1. Schematic diagram of the isulator/quantum 
dots/vacuum (IQV) tunnel junction device. 

Fig. 2. Figure of merit ZT as a function of temperature 
for various values of Ain the absence of interdot Coulomb 
interactions. Solid lines and dotted lines correspond to 
U = and U — 125T, respectively. 

Fig. 3. Figure of merit ZT as a function of temperature 
for different quantum dot configurations. 



Fig. 4. Figure of merit ZT as a function of A for 
different temperatures. 

Fig. 5. Figure of merit ZT as a function of A for 
different inelastic scattering strengthcs at ksT — 15I\ 

Fig. 6. Electrical conductance G e , thermal power S 
and electron thermal conductance Kg ctS cl function of 
temperature for different quantum dot configurations. 

Fig. 7. Electrical conductance G e , thermal power S 
and electron thermal conductance n e as a function of 
applied gate voltage for different temperatures at A = 
30T, U = 125T and U tj = 45F 

Fig. 8. Figure of merit as a function of applied gate 
voltage at k B T = 2T, A = 30T, U = 125F" and U tj = 45T. 

Fig. 9. (a) Heat current (b) average occupation num- 
ber, and (c) differential thermal conductance as a func- 
tion of AT for various values of Eb for a two-QD junc- 
tion. Y AR = 0,T] A = 0.3 and AT = 200T. 

Fig. 10. Rectification efficiency as a function of AT 
for two different values of Eb- Other parameters are the 
same as those of Fig. 9. 

Fig. 11. (a) Heat current, (b) differential thermal 
conductance and (c) thermal power as a function of AT 
for various values of Tar for a three-QD junction. 

Fig. 12. (a) Heat current, (b) differential thermal 
conductance and (c) thermal power as functions of AT 
for various values of Eq for a three-QD junction with 
T A r = and r/ A = 0.3. 

Fig. 13. (a) Heat current, (b) electrochemical po- 
tential and (c) occupation number as a function of AT 
for various values of U ac f° r a thrcc-QD junction with 
Ec = Ep + AT/5. All other parameters are same as in 
Fig. 12. 

Fig. 14. Rectification efficiency as a function of AT 
for two different values of Ec- Other parameters are the 
same as those of Fig. 12. 
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